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Введение. Исследован процесс свободных колебаний чувствительного элемента микромеханического 
гироскопа В-В-типа рамочной конструкции разработки Научно-исследовательского института прикладной 
механики им. академика В. И. Кузнецова с учетом неортогональности осей торсионов. Изучено влияние 
инструментальной погрешности изготовления на точность гироскопа на подвижном основании в случае 
свободных колебаний. Целью работы являлось повышение точности прибора посредством разработки 
математической модели микромеханического гироскопа К-К-типа с учетом неортогональности осей торсионов 
и исследования влияния указанной погрешности на точность прибора. Актуальность задачи повышения 
точности микромеханических гироскопов связана с повышением точности инерциальных навигационных 
систем, основанных на микромеханических датчиках. 

Материалы и методы. Предложены новая математическая модель, описывающая динамику гироскопа с 
учетом инструментальной погрешности изготовления прибора, и формула для оценки погрешности гироскопа. 
Приведены зависимости переменных состояния, полученные по результатам моделирования и на основе 
эксперимента. В исследовании использованы методы теоретической механики и асимптотические методы, в 
том числе формализм Лагранжа и методика осреднения Крылова-Боголюбова. 

Результаты исследования. Разработана новая математическая модель динамики гироскопа с учетом 
неортогональности осей торсионов. Получены решение уравнений малых колебаний чувствительного элемента 
гироскопа и оценка угла прецессии для случая подвижного основания. Проведен сравнительный анализ 
разработанной модели с экспериментальными данными, полученными в случае свободных колебаний 
чувствительного элемента гироскопа при неподвижном основании. По результатам анализа подтверждена 
адекватность построенной математической модели. Сформированы аналитические выражения, 
демонстрирующие тот факт, что неортогональность осей торсионов приводит к перекрестному влиянию 
амплитуд первичных колебаний на амплитуды вторичных колебаний чувствительного элемента и появлению 
дополнительной погрешности в показаниях угловой скорости при работе гироскопа в свободном режиме. 
Обсуждение и заключения. Полученные результаты могут быть использованы для повышения точности 
прибора с помощью алгоритма аналитической компенсации погрешности гироскопа и методики 
идентификации параметров математической модели. 
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Введение. Разработка высокоточных микромеханических инерциальных датчиков, в том числе 
микромеханических гироскопов (ММГ), применяемых для решения задач навигации и управления движением 
летательных аппаратов и мобильных роботов, является актуальной задачей приборостроения [1]. К 
достоинствам ММГ относятся малые масса и размеры, а также низкая стоимость по сравнению с гироскопами, 
основанными на других физических принципах. Однако к главным недостаткам ММГ относятся изменчивость 
его метрологических характеристик и низкая точность измерений параметров углового движения объекта 
(угловой скорости и угла поворота). Принцип функционирования вибрационных гироскопов основан на 
свойстве маятника Фуко сохранять плоскость малых колебаний неподвижной в инерциальном пространстве [2]. 

Основы теории гироскопов класса обобщенного маятника Фуко, к числу которых относятся и ММГ, 
изложены работах [2—5]. В них описаны различные конструктивные схемы построения ММГ, исследовано 
влияние инструментальных погрешностей изготовления и изменяющихся условий функционирования на 
динамику гироскопа. Принципиальной особенностью гироскопов класса обобщенного маятника Фуко является 
нелинейность, обусловленная конечными колебаниями чувствительных элементов (ЧЭ) или физической 
нелинейностью, связанной с особенностями системы управления колебаниями [2—6]. 

Исследования, посвященные динамике и конструированию ММГ, были опубликованы и в работах 
зарубежных авторов [6-8]. Так, например, в публикациях [6,8] получена формула для оценки уходов 
гироскопа, основанная на использовании разработанной математической модели движения, которая описывает 
медленное изменение тороидальных координат колебаний ЧЭ. В работах [0,7] обсуждаются вопросы 
изготовления ММГ и анализируются уравнения его малых колебаний. В статье [0] составляются уравнения 
движения ММГ с угловым (К-К-типа) и линейным (Г-Г-типа) осцилляторными видами движения ЧЭ. В 
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указанной статье проводится сравнительный анализ динамики таких приборов в рамках линейных моделей и 
даны рекомендации по выбору параметров ММГ, исходя из условий повышения чувствительности и 
обеспечения требуемой полосы пропускания, а также требований к линейности масштабного коэффициента. 

При проектировании ММГ разработчики стремятся использовать явление внутреннего резонанса в 
системе, обусловленное совмещением собственных частот колебаний ЧЭ [3,4]. Однако в работах [0,7] 
отмечается, что погрешности технологии изготовления, неизвестные и непредсказуемые отклонения элементов 
конструкции от проектных положений приводят к дополнительным погрешностям в измерениях прибора. 

Для повышения точности измерения угловой скорости ММГ поставлена цель: исследовать свободные 
колебания (при отсутствии управления) ЧЭ ММГ К-К-типа с учетом эффектов, возникающих из-за 
неортогональности осей торсионов. Этот дефект появляется в виду несовершенства технологии изготовления 
прибора. Ставятся задачи разработки новой математической модели динамики ММГ с учетом 
неортогональности осей торсионов, оценки ухода прибора и описания влияния неортогональности осей 
торсионов на динамику ЧЭ ММГ. 

Материалы и методы. Рассматривается модельная конструкция вибрационного ММГ В-К-типа — 
конструкция с промежуточной рамкой в соответствии с классификацией из источника [3]. Кинематическая 
схема гироскопа (рис. 1) реализована в виде двухстепенного карданного подвеса ЧЭ. 





Рис. 1. Конструктивная схема прибора: | — основание (корпус); 2 — промежуточная (внешняя) рамка; 3 — чувствительный 
элемент, состоящий из сбалансированной пластины и инерционной массы; 4 — торсионы 


Для описания положения ЧЭ введем системы координат (рис. 2), связанных: с корпусом прибора — 
ОХУХ ; с внешней рамкой упругого подвеса — Оху2,; с сбалансированной пластиной — Оху2. Причем ОЙ 
является осью чувствительности гироскопа, а система координат Ох,у,2, отличается от системы Оху2 


поворотом на постоянный угол неортогональности всех торсионов. В представленных системах начало 
координат соответствует точке Ои находится в геометрическом центре сбалансированной пластины. 
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Рис. 2. Системы координат 
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В системе Оху2, зададим осевые моменты инерции промежуточной рамки Л„, Л, Л», а в системе 
Оху; — осевые моменты инерции ЧЭ „Ло, „Л»,, Л›.. Отметим, что в данной работе оси систем координат Ох у2, и 
Оху2 считаются главными центральными осями инерции промежуточной рамки и ЧЭ соответственно. 


При моделировании движения ЧЭ принимается допущение, что конструкция торсионов обеспечивает 
бесконечную жесткость на изгиб. Положение ЧЭ относительно основания ММГ опишем двумя обобщенными 
координатами — углами а и В, а также малым постоянным углом д, характеризующим неортогональность 


осей торсионов (рис. 2). Взаимное расположение систем координат определяется последовательностью 
элементарных поворотов: 
[91 9 В 
ОХА — > ОХУ2 ——_>0%,У,2, О . 
1 2 
где под каждой стрелкой указывается ось, вокруг которой происходит поворот против хода часовой стрелки на 


угол, указанный над соответствующей стрелкой. 
Составим уравнения динамики ЧЭ ММГ в форме уравнений Лагранжа 2-го рода [9, 10]: 


4566 __0Ф 4[0\_9__0Ф 
&\0@] да да’ %\ ов) вв 9’ 


где Г, =ТГ-—// — функция Лагранжа; Ти ГП — кинетическая и потенциальная энергии системы соответственно; 


(1) 


Ф — диссипативная функция, характеризующая потери на внутреннее трение. Выражения для этих величин 
имеют вид: 
Ел + +?) + (ла? + (Лоза + Л. воз? о? 
> ›„@, + Л. „6, + Л›.6; > 1 @ „5 @-+.Л, с05°@ у 
(2) 
1..2 142 1 2 1 2 
Ф = - ОР ‚ = о” СВ ь 
где 4,, 4, — коэффициенты трения; с,, с, — коэффициенты жесткости торсионов. 
Выражения для проекций ©, ®,, в. угловой скорости ЧЭ на подвижные оси х, у, 2 имеют вид: 
@. = а с0$Всо50 — О созазтВ + зт а соз шо, 
@, =В+ 05 @ с0$0 —@зтО, (3) 


@. = азтрВсо$0 + 62 с0$ @ с0$В + ОФзш а зтВзт д. 
Учитывая малость углов @ , В и 0, тригонометрические функции в выражениях (2) и (3) от указанных 


углов можно заменить их разложениями в ряд Тейлора, ограничившись слагаемыми до первого порядка 
малости. Тогда получим уравнения малых колебаний из уравнений движения (1), записанных с учетом 
выражений (2) и (3). 

Используя формализм Лагранжа [9] в случае постоянной угловой скорости основания, получаем 
уравнения малых колебаний ЧЭ, записанные с точностью до слагаемых первого порядка малости в виде: 


.. 2 —_ с о . Л г м 2р _ | . р ^ .. 
ао. а = лОВ- —ба+— дБ, В+®,р=-лро<а--—-В+оа, (4) 
О, 2 О, 
где введены следующие обозначения (аналогично тому, как это выполнено в статье [11]: 





ВИ У. в о ЗИ >. т —4.. Со 
а а 
1х Е 2х 2у 1х т 2х 
= в О = ЗЕ: ИИ О = _@_ 
В _ 2 а _ 2 в —_ ь 
р в, (Л, +...) ЕЛ: 
Здесь 2 И — безразмерные моменты инерции упругого подвеса. , в И О О, —- собственные частоты 


колебаний и добротности по углам а , В соответственно. 

При выводе уравнений колебаний (4) угловая скорость корпуса гироскопа 6» считалась малой 
относительно собственной частоты @, , т.е. (@) << ®,, а также угол д полагался малой величиной, т. е. д <<1. 
Отметим, что в уравнениях (4) отброшены слагаемые, обусловленные наличием геометрической нелинейности 


ММГ. Влияние нелинейности геометрии движения ЧЭ на динамику ММГ В-В-типа описано в монографии [5]. 
Принимая во внимание, что правые части уравнений (4) являются малыми возмущениями, т. е. 


5 2 ;. 2 
й+®,@=О(=), то с точностью до слагаемых первого порядка малости можно записать: & =-®.@+0О(=). 


Таким образом, из правой части уравнений (4) исключаются вторые производные углов @ иВ. 
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Рассмотрим случай изотропного упругого подвеса, т.е. равенства собственных частот колебаний и 
равных добротностей: 


@„, — @, = @, О, =О, =0, 
где &, — характерное значение собственной частоты колебаний; () — характерное значение добротности. 
Отметим, что случай наличия разнодобротнсти (0, 5 О,) и малой разночастотности (®, ге а, при 


исследовании свободных колебаний ЧЭ ММГ рассмотрен в статье [11]. При введенных обозначениях и 
принятых допущениях запишем уравнения движения ЧЭ с точностью до слагаемых первого порядка малости в 
безразмерном виде: 


Я-+ о = ЛОВ _ О‘ а-= —— т ев, 
Л (5) 
В+о2В =-рОЯ-О"'®,В- бога. 

Отметим, что система уравнений (5) приводится к стандартной форме записи регулярно возмущенной 
системы дифференциальных уравнений с одной быстрой угловой переменной [12, 14]. Одним из 
распространенных способов нахождения решения регулярно возмущенных систем является использование 
асимптотических методов разделения движений [12—15]. 

Решение нелинейных уравнений (5) получим с помощью методики осреднения Крылова — 
Боголюбова [13], а в качестве медленно меняющихся переменных будем использовать переменные Ван-дер- 
Поля [12] р, Чл, Р>, 42: 

а = рзш(®,Ё)-+ 4, с03(61), —@а=вур, с03(®1)-®, 4, зт(®!), 
В= р, зщ (0/) +4, с0$(®1), — В=ор, соз (в) -®,а, зт (1). 


С помощью процедуры осреднения [14,15] по явно входящему времени, получаем осредненную 
систему дифференциальных уравнений, разрешенную относительно производных медленных переменных Ван- 











дер-Поля: 
| | р) 7) | ло И) 
р=-=-О"р+я 4. а=--О 4+5 += р., 
— 24, 2.) 2 @) А (6) 
р 1 —1 РО 0 р 1 —1 э 9 
НЫ — ей, Зы — +—р. 
р> 59 р> рум Р! 54 (> 59 (> 5 5 р 


Штрихом в уравнениях (6) обозначено дифференцирование по безразмерному времени т =@. 


Полученная модель в виде линейных дифференциальных уравнений описывает свободные колебания 
ЧЭ гироскопа на подвижном основании. Решение системы уравнений (6) можно записать в виде: 


р (т) =ехр в Росоз[ № +7?) дров) [| +7 "г , 
Ур Му +7 


4, (т) = ехр 5 в соз( и? +7?) + д Рут , 


(7) 
р» (т) =ехр 5 Рх соз( У’ +у МИ +=) ее ы( фир | 
20 и + 
4, (т) = ехр 5 Ч сз У’ +У Ми +иг)- ра о Ры) и +7? =) 
20 +? 
гле рь=р (0), р» =р, (0), а, =а4, (0), 4» =а, (0) — начальные условия; безразмерная угловая скорость 


основания прибора и=.| ий »0/(2%,); параметр, характеризующий неортогональность осей торсионов 


7-52 


Вторые слагаемые в формулах (7) характеризуют перекрестное влияние первичных колебаний на 
вторичные колебания и наоборот. Отметим, что в случае ортогональных осей торсионов, когда у=0, 


решение (7) совпадает с результатами работы [5]. 

Полученные аналитические решения (7) уравнений колебаний представляют интерес для разработки 
методик идентификации параметров, а также прогнозирования ухода гироскопа и его учета при использовании 
методики алгоритмической компенсации погрешностей. 
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Результаты исследования. Для валидации разработанной модели сравним результаты моделирования, 
вычисленные по формулам (7), с экспериментальными данными. Измерительная информация получена с 
помощью системы наблюдения. В качестве измерительной информации электростатических датчиков имеем 
переменные Ван-дер-Поля ру, 41, ро, 42. 

В эксперименте использован образец прибора со следующими параметрами математической модели 


при неподвижном основании (у=0):О=3856, л=),=ЪЬ у=0,2.10”, в качестве начальных условий для 
переменных Ван-дер-Поля выбраны значения равные измерениям в начальный момент времени: 
Ро =13,467-10°, а, =20,429.10°, р» =0,787-10”, а» =1172.10>. 
Значение параметра у соответствует углу неортогональности осей торсионов д равного одной 


угловой секунде. Графическое представление зависимостей переменных Ван-дер-Поля р,(т), а, (т), р. (т), 


4, (г) медленно меняющихся с течением безразмерного времени, приведено на рисунке 3. 


р1-10? р2-10° 
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| 
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Рис. 3. Зависимости переменных Ван-дер-Поля: | — результаты моделирования; 2 — экспериментальные данные 


Графики (рис. 3) демонстрируют значительное совпадение зависимостей для переменных р,, 4,, 


полученных По результатам моделирования, с экспериментальными данными. Зависимости для переменных 


Р. Ч. Полученные по результатам моделирования, качественно согласуются с данными эксперимента, а 


наблюдаемые небольшие количественные отклонения могут быть обусловлены нелинейными эффектами, 
такими как нелинейность геометрии движения ЧЭ [5], или явлениями разнодобротности, разночастотности и 
погрешностью смещения инерционной массы [11]. Рассмотрение нелинейных эффектов, влияющих на 
динамику ММГ, при построении математических моделей колебаний ЧЭ увеличивает точность 
микромеханических датчиков в составе инерциальных навигационных систем [2]. 

Уход гироскопа из-за нелинейных эффектов и других инструментальных погрешностей будем 
оценивать с помощью вспомогательного функционала / [5, 6, 8]: 


1 2УЛЛ (94 +РР.) | 
л(а+р)-л(е-+в) 


который связан с углом 9 через соотношение: 


(8) 


9= 5хагаал (1) | 


Причем этот параметр пропорционален интегралу от угловой скорости: 
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С помощью формулы (8), учитывая решение (7), можно оценивать уход гироскопа, связанный с 


9 = ЕЕ Го(=)аз | 
0 


неортогональностью осей торсионов, возникшей из-за несовершенства технологии изготовления. На рис. 4. 
представлены зависимости функционала / от безразмерного времени по результатам эксперимента и расчета 


по формуле (8). 
[ 
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Рис. 4. Зависимость функционала / (т) : | — результаты моделирования; 2 — экспериментальные данные 


Близость оценки ухода (рис. 4) и переменных Ван-дер-Поля (рис. 3), полученных аналитически, с 
экспериментальными данными характеризуют хорошую точность построенной модели, особенно если учесть 
тот факт, что в ней пренебрегались нелинейные эффекты, а также явления разночастотности и 
разнодобротности. Несмотря на указанные пренебрежения, модель позволяет построить методики 
идентификации параметров, с помощью которых можно уточнить зависимости, полученные при 
моделировании. Применение методик идентификации параметров математической модели приведет к 
повышению точности ММГ в режиме вынужденных колебаний, который является рабочим режимом 
гироскопов. 

Обсуждение и заключения. Построена новая математическая модель ММГ К-К-типа для режима 
свободных колебаний ЧЭ. В модели учитывается неортогональность осей торсионов, возникающая в следствие 
технологической невозможности обеспечить высокую точность изготовления прибора. Получена формула 
оценки угла прецессии при подвижном основании прибора. Посредством сравнения результатов моделирования 
с данными эксперимента проведена валидация математической модели ММГ. Показано, что неортогональность 
осей торсионов приводит к перекрестному влиянию первичных колебаний на величину вторичных колебаний и 
наоборот. Результаты работы могут быть использованы в алгоритме аналитической компенсации погрешности 
гироскопа с целью повышения точности ММГ. 
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